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Abstract 

We explore electron transport properties in molecular wires made of heterocyclic molecules (pyrrole, furan 
and thiophene) by using the Green's function technique. Parametric calculations are given based on the 
tight-binding model to describe the electron transport in these wires. It is observed that the transport 
properties are significantly influenced by (a) the hctcroatoms in the heterocyclic molecules and (b) the 
molecule-to-electrodes coupling strength. Conductance (g) shows sharp resonance peaks associated with 
the molecular energy levels in the limit of weak molecular coupling, while they get broadened in the 
strong molecular coupling limit. These resonances get shifted with the change of the heteroatoms in 
these heterocyclic molecules. All the essential features of the electron transfer through these molecular 
wires become much more clearly visible from the study of our current-voltage (I-V) characteristics, and 
they provide several key informations in the study of molecular transport. 
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1 Introduction 

Quantum transport properties of organic molecules 
bridging over electrodes are recent interest of nan- 
otechnologies since they constitute promising build- 
ing blocks for future generation of nanoclcctronic 
devices where the electron transport is predomi- 
nantly coherent. [U [2] Following experimental devel- 
opments, theory can play a major role in under- 
standing the new mechanisms of conductance, and 
in the last few decades, electron transport through 
different nano-scale systems [31 HI G3 [5] have been 
studied enormously. The single-molecule electron- 
ics plays a key role in the design of future nano- 
electronic circuits, but the goal of developing a re- 
liable molecular-electronics technology is still over 
the horizon and many key problems, such as de- 
vice stability, reproducibility and the control of 
single-molecule transport need to be solved. Start- 
ing with the paper of Aviram and Ratner[7] in 
which a molecular electronic device has been sug- 
gested for the first time, the development of a 
theoretical description of molecular electronic de- 
vices has been pursued. Since then several nu- 
merous experiments [9j QUI [UJ [12] have been 
performed through molecules placed between two 
metallic electrodes with few nanometer separation. 
The operation of such two-terminal devices is due 
to an applied bias. Current passing across the 
junction is strongly nonlinear function of the ap- 
plied bias voltage and its detailed description is 
a very complex problem. The complete knowl- 
edge of the conduction mechanism in this scale 
is not well understood even today. For exam- 
ple, it is not very clear how the molecular trans- 
port is affected by the structure of the molecule 
itself or by the nature of its coupling with the 
electrodes. The most important issue is proba- 
bly the quantum interference effects [131 [TH [T51 [To] 
among the electron waves traversing through dif- 
ferent arms of the molecule. Another important 
issue is the molecular coupling to the side attached 
electrodes. [17] Tuning this coupling, one can control 
the current amplitude across the bridge quite sig- 
nificantly. Similar to these, there are several other 
factors like the electron-electron correlations. p~8] 
dynamical fluctuations. [191 l20j etc., which provide 
rich effects in the electron transport. To design 
molecular electronic devices with specific proper- 
ties, structure-conductance relationships are also 
needed, and in a very recent work Ernzerhof et 
aL[2"T] have presented a general design principle 
and performed several model calculations to demon- 



strate the concept. Here we focus on single molecu- 
lar transport that are currently the subject of sub- 
stantial experimental, theoretical and technologi- 
cal interest. These molecular systems can act as 
gates, switches, or transport elements, providing 
new molecular functions that need to be well char- 
acterized and understood. In many molecular de- 
vices, electron transport is dominated by conduc- 
tion through broadened HOMO or LUMO states. 
In contrast, in this article we find that the elec- 
tron transport through molecular bridges can be 
controlled very sensitively by chemically modifying 
the heteroatoms in the heterocyclic molecules. 

Here we demonstrate that there are advantages 
in engineering molecules which can be modified ex- 
ternally to achieve control over transport by alter- 
ing the properties of the heteroatoms. In an effi- 
cient molecular transport system, actual contact of 
the molecule to both electrodes is required. There- 
fore the simplest theoretical view is based on the 
tight-binding type one-electron picture. In this ar- 
ticle we reproduce an analytic approach based on 
the tight-binding model to investigate the electron 
transport properties for the model of single hetero- 
cyclic molecules named as pyrrole, furan and thio- 
phene, and the coupling of the molecules to the side 
attached electrodes is treated through the Newns- 
Anderson chemisorption theory. [22"1 [23[ [24] There 
exist several ah initio methods for the calculation 
of the conductance [23 HQ |571 HH H21 SOI EH EH 
[31 ES] as well as model calculations. [23 ] [24 l l36 l [37] 
The model calculations are motivated by the fact 
that the ah initio theories are computationally too 
expensive, while the model calculations by using the 
tight-binding formulation are computationally very 
cheap and also provide a worth insight to the prob- 
lem. In our present study, attention is drawn on 
the qualitative behavior of the physical quantities 
rather than the quantitative ones. Here we do not 
take into account the effect of charge transfer from 
the electrode which may play a significant role in 
the ah initio study. 

Our scheme of study is as follows. In Section 
2 we describe very briefly about the methodology 
for the calculation of the transmission probability 
(T) and the current (/) through a finite size con- 
ducting system sandwiched between two metallic 
electrodes by the use of Green's function technique. 
Section 3 investigates the behavior of the conduc- 
tance (g) as a function of the injecting electron en- 
ergy (E) and the current-voltage (J-V) character- 
istics for the model of three different heterocyclic 
molecules. Here we focus our results in the aspects 
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of (a) the heteroatoms in the heterocyclic molecules, 
and (b) the molecular coupling to the side attached 
electrodes. Finally, we summarize our results in 
Section 4. 



2 A Brief Description Onto 
the Theoretical Formulation 

In this section we discuss very briefly about the 
methodology for the calculation of the transmis- 
sion probability (T), conductance (g) and current 
(/) through a finite size conductor attached to two 
one-dimensional semi-infinite metallic electrodes by 
the use of Green's function technique. 

Let us refer to Fig. [TJ where a one-dimensional 
conductor with N number of atomic sites (array of 
filled circles) connected to two semi-infinite metallic 
electrodes, viz, source and drain. The conducting 
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Figure 1 : Schematic view of a one-dimensional con- 
ductor with N number of atomic sites (filled cir- 
cles) attached to two metallic electrodes (source and 
drain) through the atomic sites 1 and N respec- 
tively. 

system within the two electrodes can be anything 
like a single molecule, or an array of few molecules, 
or an array of some quantum dots, etc. At suffi- 
cient low temperature and applied bias voltage, the 
conductance of the conductor can be written by the 
Landaucr conductance formula|38, 391 as, 
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where g is the conductance and T is the transmis- 
sion probability of an electron through the conduc- 
tor. The transmission probability can be expressed 
in terms of the Green's function of the conductor 
and its coupling to the side attached electrodes by 
the relation. 1551 1591 



T = Tr [T s G r c T D G«] 



(2) 



where G r c and G° are respectively the retarded and 
advanced Green's function of the conductor. T$ 
and Yd are the coupling terms due to the coupling 
of the conductor to the source and drain respec- 
tively. For the complete system i.e., the conductor 



with the two electrodes, the Green's function is de- 
fined as, 

G=(e-H)~ 1 (3) 

where e = E + ir). E is the energy of the source 
electron and r\ gives an infinitesimal imaginary part 
to e. Evaluation of this Green's function requires 
the inversion of an infinite matrix as the system 
consists of the finite conductor and the two semi- 
infinite electrodes. However, the entire system can 
be partitioned into sub-matrices corresponding to 
the individual sub-systems, and the Green's func- 
tion for the conductor can be effectively written as, 



G c = (e - H c 



(4) 



where H c is the Hamiltonian of the conductor sand- 
wiched between the two electrodes. The Hamilto- 
nian of the conductor in the tight-binding frame- 
work can be written within the non-interacting pic- 
ture in this form, 
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where c\ (cj) is the creation (annihilation) opera- 
tor of an electron at site i, e^'s are the site energies 
and t is the nearest-neighbor hopping strength. In 
Eq. (gD, £5 = h\ c g s h Sc and E D = h Dc g D h Dc are 
the self-energy operators due to the two electrodes, 
where gs and go are respectively the Green's func- 
tion for the source and drain. hs c and hoc are the 
coupling matrices and they will be non-zero only 
for the adjacent points of the conductor, 1 and N 
as shown in Fig. [TJ and the electrodes respectively. 
The coupling terms Y$ and Yd for the conductor 
can be calculated through the expression, 
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where ^gm) ana - ^S(D) are the retarded and ad- 
vanced self-energies respectively and they are con- 
jugate with each other. Datta et al. [40] have shown 
that the self-energies can be expressed like as, 
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where Ag(m are the real parts of the self-energies 
which correspond to the shift of the energy eigen- 
values of the conductor, and the imaginary parts 
Ag^D) °f the self-energies represent the broaden- 
ing of these energy levels. This broadening is much 
larger than the thermal broadening, and accord- 
ingly, we restrict our all calculations only at abso- 
lute zero temperature. Here we adopt the Newns- 
Anderson chemisorption model|22l (23j El] for the 
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description of the electrodes and for the interaction 
of the electrodes with the conductor, where the ef- 
fect of the electrodes is then formally incorporated 
through the self-energies Eg and T,d- To describe 
the electrodes (in the form of a semi-infinite one- 
dimensional chain), in this present scheme, we use 
the similar kind of tight-binding model Hamilto- 
nian as presented in Eq. l[5]). which is parametrized 
by the constant on-site potential eg an d nearest- 
neighbor hopping integral v. From the standpoint 
of the band theory or molecular orbital, the cou- 
pling between the conductor and the electrodes can 
be attributed to one-electron hopping processes, 
where the hopping parameters are rg and tjj respec- 
tively. All these are the essential parameters for this 
particular scheme to describe the electron transport 
phenomena in a bridge system. Now by utilizing 
the Newns-Anderson type model, we can express 
the conductance in terms of the effective conductor 
properties multiplied by the effective state densities 
involving the coupling. This permits us to study 
directly the conductance as a function of the prop- 
erties of the electronic structure of the conductor 
within the bridge. 

Since the coupling matrices hs c and h]j c are non- 
zero only for the adjacent points in the conductor, 
1 and N as shown in Fig. [TJ the transmission prob- 
ability becomes, 

T(E, V) = 4Af 1 (S, V)A% N (E, V)\G 1N (E, V)\ 2 

(8) 

where An =< 1|A|1 >, A NN =< N\A\N > and 
Gin =< 1|G C |JV>. 

The current passing across the conductor is de- 
picted as a single-electron scattering process be- 
tween the two reservoirs of charge carriers. The 
current-voltage relation is evaluated from the fol- 
lowing expression. [38] 

fE F +eV/2 

I{V) = — T{E,V)dE (9) 

™ JE F -eV/2 

where Ep is the equilibrium Fermi energy. For the 
sake of simplicity, here we assume that the entire 
voltage is dropped across the conductor-electrode 
interfaces and this assumption doesn't greatly af- 
fect the qualitative aspects of the I-V character- 
istics. This assumption is based on the fact that 
the electric field inside the conductor, especially for 
short conductors, seems to have a minimal effect 
on the conductance-voltage characteristics. On the 
other hand for quite larger conductors and high bias 
voltages, the electric field inside the conductor may 
play a more significant role depending on the in- 
ternal structure and size of the conductor. [40] but 



yet the effect is too small. Using the expression of 
T(E, V) as in Eq. © the final form of I(V) will be, 



I(V) = 



4e 



E F +eV/2 
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x\G 1N (E,V)\ 2 dE 
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Eq. (P), Eq. © and Eq. ((TDJ) are the final work- 
ing formule for the calculation of the conductance 
g, transmission probability T, and current I respec- 
tively through any finite size conductor sandwiched 
between two metallic electrodes. 

In this paper, we will describe the electron trans- 
port properties by using the above methodology for 
the three different models of heterocyclic molecules 
those are defined as pyrrole, furan and thiophene 
(Fig. [2|). For simplicity, we take the unit c — e = 
h = 1 in our present calculations. 

3 Results and their Interpre- 
tation 

This section investigates the behavior of the con- 
ductance g as a function of the injecting elec- 
tron energy E, and the variation of the current / 
with the applied bias voltage V for the three dif- 
ferent molecular wires containing the heterocyclic 
molecules. Here we focus our results on the electron 
transport properties considering the effects of (a) 
the hcteroatoms in the heterocyclic molecules (pyr- 
role, furan and thiophene), and (b) the molccule- 
to-electrodes coupling strength. The arrangements 
of the three different molecular wires are shown in 
Fig. [21 where nitrogen, oxygen and sulfur are the 
heteroatoms in the molecular wires consisting with 
pyrrole, furan and thiophene molecules respectively. 
These molecules are attached to the semi-infinite 
metallic electrodes by thiol (sulfur-hydrogen i.e., S- 
H bond) groups. In actual experimental set-up, the 
electrodes made from gold (Au) are used and the 
molecule coupled to the electrodes through thiol 
groups in the chemisorption technique where hydro- 
gen (H) atoms remove and sulfur (S) atoms reside. 
To describe these heterocyclic molecules, we use the 
similar kind of non-interacting tight-binding Hamil- 
tonian as given in Eq. ([5]). 

All the essential features of this article are stud- 
ied in the two distinct regimes. One is so-called 
the weak-coupling regime, defined by the condi- 
tion T{s,d} « t. The other one is so-called the 
strong-coupling regime, denoted by the condition 
t {S,d} ~ t. For these two limiting cases, we take 
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energy values, while they (resonance peaks) disap- 
pear for all other energies. These resonance peaks 
are associated with the energy eigenvalues of the 
individual molecules. Therefore we can say that, 
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Figure 2: Schematic view of the three different 
molecular wires. The heterocyclic molecules (pyr- 
role, furan and thiophene) are attached to the elec- 
trodes, viz, source and the drain via thiol (S-H) 
groups. 



the values of the different parameters as follows: 
ts = to = 0.5; t = 2.5 (weak coupling) and 
ts = Td = 2; t = 2.5 (strong-coupling). Here we 
set the on-site energy eo = (we can take any con- 
stant value of it instead of zero since it gives only 
the reference energy level) for the electrodes, and 
the hopping strength v — 4 in the two semi-infinite 
metallic electrodes. For the sake of simplicity, here 
we set the Fermi energy Ep = 0. 

In Fig. [21 we plot the conductance g as a func- 
tion of the injecting electron energy E for the three 
molecular wires, where the solid, dotted and dashed 
curves correspond to the results for the wires con- 
taining pyrrole, furan and thiophene molecules re- 
spectively. The results for the weak-coupling case 
are shown in Fig. EJa), while Fig. 0(b) gives the re- 
sults in the limit of strong molecular coupling. In 
the weak molecular coupling, the g-E characteris- 
tics show sharp resonance peaks for some particular 



Figure 3: Conductance g as a function of the inject- 
ing electron energy E, where the solid, dotted and 
dashed curves correspond to the molecular wires 
with pyrrole, furan and thiophene molecules re- 
spectively, (a) weak-coupling limit and (b) strong- 
coupling limit. 

the conductance spectrum manifests itself the elec- 
tronic structure of the molecule. As expected, the 
maximum value of these conductance peaks goes to 
two i.e., the transmission probability (T) becomes 
unity since we get the relation g = 2T from the Lan- 
dauer conductance formula, Eq. ([1]) , with e = h = 1 
in our present formulation. From these results it is 
observed that the resonance peaks get shifted in the 
scale of the energy E as we chemically change the 
heteroatoms in these heterocyclic molecules. Thus 
we can predict that the on/off state of the elec- 
tron conduction across the molecule can be tuned 
by chemically modifying the heteroatoms, though 
in all these three molecular wires the molecular 
structures are the same. This provides an impor- 
tant finding in the study of molecular transport. 
In the strong molecule-to-electrodes coupling limit, 
these resonance peaks get substantial widths as pre- 
sented in Fig'.^b). Such increment of the resonance 
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widths is due to the broadening of the molecular en- 
ergy levels in the limit of strong molecular coupling, 
where the contribution comes from the imaginary 
parts of the self-energies Eg and Ed, [38] as men- 
tioned earlier in Section 2. From the curves plotted 
in Figs. [3ja) and (b) it is observed that, the po- 
sitions of the resonance peaks are independent of 
the molecule-to-electrodes coupling strength. An- 
other significant feature observed from these curves 
is that, in the strong-coupling limit, the molecular 
bridge remains in the on-state condition i.e., elec- 
tron passes through the molecule for a wide range 
of energies, while a fine tuning in the energy scale 
is necessary to get the on-state condition for that 
bridge in the limit of weak molecular coupling. This 
feature is quite significant in fabrication of efficient 
electronic circuits by using these molecules. 

The scenario of the electron transfer through such 
molecular wires becomes much more clearly visi- 
ble by describing the current-voltage characteris- 
tics, where the current passing across the molecule 
is computed from the integration procedure of the 
transmission function T (see Eq. ©). The variation 
of this transmission function looks similar to that 
of the conductance spectra (see Figs.EJa) and (b)), 
differ only in magnitude by the factor 2, since the 
relation g = 2T holds from the Landauer conduc- 
tance formula as stated earlier. In Fig.|H we display 
the current I as a function of the applied bias volt- 
age V for the three molecular wires where the solid, 
dotted and dashed curves correspond to the same 
molecular systems as given in Fig. |3l The results 
for the weak-coupling limit are shown in Fig. |U[a) , 
while Fig. 0Jb) corresponds to the results in the 
limit of strong molecular coupling. In the limit of 
weak molecular coupling, the current-voltage char- 
acteristics give staircase-like behavior with sharp 
steps (Fig. (Ha)). This is due to the existence of 
the sharp resonance peaks in the conductance spec- 
trum (Fig. [3ja)) since the current is evaluated from 
the integration procedure of the transmission func- 
tion T. With the increase of the applied bias volt- 
age, the electrochemical potentials on the electrodes 
are shifted gradually, and eventually cross one of 
the molecular energy levels. Accordingly, a current 
channel is opened up and the current-voltage curve 
produces a jump. The shape and height of these 
current steps depend on the width of the molecular 
resonances. With the increase of the molecular cou- 
pling strength, the current varies quite continuously 
with the bias voltage V, as illustrated in Fig. Sfb), 
and achieves much higher values. This can be un- 
derstood by noting the areas under the curves in 



the conductance spectrum as plotted in Fig. Efb) 
in the limit of strong molecular coupling. From 
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Figure 4: Current J as a function of the applied 
bias voltage V, where the solid, dotted and dashed 
curves correspond to the molecular wires with pyr- 
role, furan and thiophene molecules respectively, 
(a) weak-coupling limit and (b) strong-coupling 
limit. 

these results we predict that one can achieve much 
higher current amplitude in a molecular bridge sys- 
tem just by tuning the molecule-to-electrodes cou- 
pling strength, without changing any geometry of 
the bridge. The other most important feature ap- 
pears from these current-voltage characteristics is 
that, the threshold bias voltage (Vth) of the elec- 
tron conduction changes with the change of the 
heteroatoms in these heterocyclic molecules. Thus 
we can tune Vth by chemically modifying the het- 
eroatoms and this provides a key result in the study 
of molecular transport in these molecular wires. 

4 Conclusion 

In conclusion, we have explored the electron trans- 
port properties of the molecular wires made by 
the heterocyclic molecules (pyrrole, furan and thio- 
phene), and observed that the transport properties 
are significantly influenced by (a) the heteroatoms 
in the molecules and (b) the molecule-to-electrodes 



6 



coupling strength. Here we have used the simple 
tight-binding model to describe the molecular wires, 
and introduced a parametric approach to study the 
electron transport. 

In the study of the g-E characteristics we have 
noticed that the conductance shows very sharp res- 
onance peaks (Fig.[3^a)) in the limit of weak molec- 
ular coupling associated with the molecular en- 
ergy levels, while the widths of these resonances 
get enhanced substantially (Fig.Efb)) in the strong 
molecular coupling limit. This is due to the broad- 
ening of the molecular energy levels in the limit of 
strong coupling, where the contribution comes from 
the imaginary parts of the two self energies Eg and 

Next we have concentrated our study on the 
current-voltage characteristics from which the sce- 
nario of the electron transfer through the molecules 
can be understood much more clearly. The current 
changes its behavior from the staircase-like struc- 
ture with sharp steps (Fig. |4ja)) to the continuous 
one (Fig. 0Jb)) as the molecular coupling changes 
its strength from the weak regime to the strong one. 
From this study we have observed that the thresh- 
old bias voltage Vth of the electron conduction can 
be tuned by chemically modifying the heteroatoms 
in these heterocyclic molecules. 

Throughout our study we have used several re- 
alistic assumptions. More studies are expected to 
take the Schottky effect, comes from the charge 
transfer across the metal-molecule interfaces, the 
static Stark effect, which is taken into account for 
the modification of the electronic structure of the 
molecular bridge due to the applied bias voltage 
(essential especially for higher voltages). However 
all these effects can be included into our framework 
by a simple generalization of the given formalism 
described here. In this article we have also ignored 
the effects of all inelastic scattering processes, by as- 
suming that the electrons move smoothly from the 
source to the drain subject only to elastic scattering 
within the junction, and electron-electron correla- 
tion to characterize the electronic transport through 
such heterocyclic molecules. 
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